Scientific Reports

Central to understanding biodiversity is the systematic classification of biological taxa where concepts on integrating omics data have not yet been comprehensively executed. The data-driven approach is aiming at integrating phylogenetic data from three or more scales aiming at enriching the robustness of species delimitations while also revealing evolutionary drivers of diversification and speciation. This study investigates complex-thallose liverworts as reference and combines the analysis of DNA marker sequencing, morphometrics utilizing bioimaging measurements and chemometrics using liquid chromatography high-resolution mass-spectrometry (UPLC/ESI-QTOF-MS) with data-dependent acquisition of tandem mass-spectra (DDA-MS), with character variation evaluated with dendrograms and data mining. The comparative analysis revealed similar tree topologies albeit with diverging positions of taxa that were attributed to flavonoid-glycosides, auronidins and phenanthrenes resulting from evolutionary radiation and adaptations to stressful bioclimatic conditions, whereas unique fatty acyls in Riccia link to ecological causes of speciation. Characteristic amino acid motifs confirm the intermediary position of liverworts between algae and land plants. By functionally attributing and mapping chemometric markers to the taxonomy, the data-driven approach integrating omics into integrative taxonomy is opening opportunities for plant systematics to identify evolutionary drivers and form new hypotheses on the diversification and speciation of species.

Age and sex estimation from dental radiographs supports forensic identification and age-dependent legal assessment. Multi-input models may combine global and regional representations, but when these inputs are derived from the same radiograph, they may repeat anatomical information rather than provide complementary signals. We developed HMA-Net, a hierarchical multi-stream anatomical network, together with a task-specific input audit that evaluates whether each candidate representation improves at least one prediction task without materially degrading the other. The model was developed using a retrospective single-centre cohort of 4,293 panoramic radiographs from individuals aged 3-18 years and combined a panoramic-image stream with graph-linked tooth representations. On the held-out test set, the final panoramic-plus-permanent-tooth model achieved a mean absolute age error of 0.51 years and a sex-classification accuracy of 85.2%; 99.4% of age estimates were within 2 years of chronological age. Across three independent training runs, removing the jaw-region stream reduced mean absolute age error from 0.68 +/- 0.01 to 0.52 +/- 0.01 years and increased sex accuracy from 81.3+/- 1.7% to 85.8 +/- 0.8%. Adding 20 deciduous-tooth classes provided no age benefit and reduced sex accuracy by 7.9% relative to the permanent-only representation. The mandibular first molars received the greatest tooth-level attention, and their occlusion impaired both prediction tasks. These findings show that anatomically plausible inputs do not necessarily provide incremental predictive value. Task-specific input auditing can improve joint model performance, reduce input redundancy and preserve anatomically traceable predictions. More broadly, this framework provides a testable approach for selecting correlated inputs in medical-imaging systems in which regional views re-present anatomy already contained in the whole image. Keywords: Forensic odontology; Dental age estimation; Sex classification; Panoramic radiography; Multi-input learning; Anatomical input auditing

The erythrocyte sedimentation rate (ESR) is one of the most common and widely used laboratory diagnostic parameters in connection with inflammatory reactions and it is probable that every reader has already experienced a determination of their ESR. A rapid ESR is a non-specific parameter that provides information about the inflammatory process. Although the origins of this methodology date back to antiquity, the description of the process as the collapse of a percolating gel formed from erythrocytes has only recently been achieved. It was not yet known whether slow ESR has any medically relevant significance. Here we show a variety of clinical pictures that exhibit a systematically slow ESR (e.g., sickle cell disease, neuroacanthocytosis syndromes, chronic mountain sickness). Using a combination of measured data and physical modelling, we show how the accuracy and significance of ESR data can be increased. With this improved ESR (supraESR), we introduce a completely new, cost-effective diagnostic parameter, based on an established and easily automated measurement method, that enables low-cost screening for neuroacanthocytosis syndrome, a group of rare neurodegenerative diseases previously detectable only through complex diagnostic tests.

Single cell differential expression analysis enables biologists to make statistical conclusions about which genes are up or downregulated in a particular cell type, between two conditions, such as those with or without a disease. However, due to biological and technical noise, these differences are hard to detect reliably. Determining if an experiment has sufficient power is often derived from simulations or small pilot samples, if done at all. Here, we use sex-biased differential expression on 1,494 donors in three brain cell types to derive empirically grounded power estimates for a variety of experimental setups. Our work reveals a substantial lack of power for reliably detecting the small effects in the range that many studies report, even in experiments containing 600 donors. When reducing the astrocyte data to the poor sequencing characteristics of microglia, over half the differentially expressed genes (DEGs) detected in the full set were lost, highlighting the damage caused by insufficient sequencing depth and cell counts. Despite standard thresholds of adjusted p-value with multiple testing correction, only the top quartile of significant genes were reproducible. Furthermore, from fitting a predictive model to the empirical outcomes, we find cell count and expression levels of genes to be a strong determinant of power. As such, we advocate for future studies to employ cell type enrichment and deeper sequencing, especially for rarer populations like microglia, emphasise the importance of powering experiments of this type when seeking robust findings, and suggest stricter significance thresholds for future discoveries.

Dominance hierarchies are a fundamental aspect of social organization, enabling animals to minimize aggression and optimize access to resources. Previous studies have highlighted the energetic and physiological demands of dominant status, as well as the behavioral flexibility required of subordinates to navigate these hierarchies. Despite advancements in automated behavior tracking, limitations persist in tracking fine-scale, real-time interactions within complex social environments. Here, we developed and validated a novel RFID-based system to continuously monitor dominance hierarchies in group-housed male mice over 10 days. This system enabled unbiased behavioral inference across light phases and revealed spatial and temporal patterns of dominance behavior undetectable through traditional live-scored methods. Automated tracking accurately identified alpha individuals and consistently inferred linear hierarchies across cohorts, with greater precision for higher-ranked individuals. Behavioral metrics, such as transition frequencies and proximity to food zones, were consistent with dominance driven activity. Hormonal analyses revealed that higher-ranked mice exhibited increased leptin and peptide YY, consistent with heightened activity and satiety signaling, while lower C-peptide levels reflected greater metabolic demands of dominance. Furthermore, dominance rank was associated with differences in light-dark activity, which were in turn related to circulating hormone profiles. This study demonstrates the utility of automated RFID tracking in capturing dominance hierarchies with temporal and spatial granularity, while revealing links between social rank, metabolic regulation, and activity patterns advancing our understanding of social behavior dynamics.

Early-life social experiences with parents and littermates are critical for the development of social behavior and cognition in social species. Yet, how naturalistic social interactions evolve across development remains poorly understood, because reliable tracking of multiple freely interacting animals is challenged by frequent occlusion and body overlap. To address this challenge, we developed Move-Altogether (MovAl), a multi-animal pose estimation pipeline integrating video object segmentation and contour-guided preprocessing that improved identity stability and pose accuracy compared with widely used methods. Using MovAl-derived pose sequences, we established a denoising autoencoder (DAE)-based framework to map longitudinal home-cage behavior within a common three-dimensional space defined by inter-animal spacing, social orientation, and interaction dynamics. Within this space, control mice showed progressive reorganization of group behavior from 3 to 8 weeks of age, establishing a developmental reference for naturalistic social organization. Against this reference, prenatal valproic acid (VPA) exposure produced prominent alterations in social behavior organization at 3-5 weeks of age, demonstrating that the consequences of neurodevelopmental perturbation can be detected through changes in naturalistic group behavior among littermates. Extending the analysis to preweaning mother-pup interactions revealed increased maternal movement and reduced mother-pup proximity in the VPA condition, indicating alterations in the early social environment that preceded postweaning behavioral differences. These findings suggest that quantitative characterization of early social interactions may provide a means to identify early behavioral alterations associated with atypical neurodevelopment. Together, the MovAl-DAE framework enables longitudinal quantification of naturalistic social behavior, defining its developmental organization and revealing alterations associated with neurodevelopmental perturbation.

Kidney diseases are a significant health concern worldwide, with a growing demand for fundamental treatments, such as generation of transplantable organs from programmed cells. The generation of kidney organoids has led to significant advances in disease modeling. However, failure to establish a robust endothelial vasculature network in kidney organoids remains a barrier to oxygen support and gas exchange in renal tissue engineering. We targeted the development of early kidney vasculature system by advancing the classic mouse embryonic kidney organ culture and organoid model system by introducing a novel mini reservoir setup. The method improved integration with the highly vascularized chicken Chorion Allantois Membrane (CAM), enabling host blood flow to developing kidney organoids. We examined the origin and role of the CD31+ endothelial cells (ECs) and their CD146+ progenitors by applying a dissociated embryonic kidney approach. This enabled selected depletion of specific cell populations, followed by regeneration of the kidney organ primordia by reaggregation and placement on the CAM to assess organogenesis. The cell lineage and fate tracing with constitutively labelled cells revealed that the intrinsic endothelial CD146+ progenitor and those of CD31+ ECs in the kidney initiating organogenesis established the kidney vasculature. If the intrinsic EC progenitors were depleted, the mouse glomerular vascularization was rescued by the CAM derived ECs.

Background: Primary embryonic dorsal root ganglion (DRG) explants offer an accessible and physiologically relevant organotypic in vitro assay platform that preserves native three-dimensional cellular interactions and supports robust neurite outgrowth, enabling studies of sensory neuron development and neurite outgrowth, as well as controlled assessment of neurotoxicity, and pharmacological compound effects. Methods: Here we provide a comprehensive step-by-step method to micro-dissect and culture embryonic chicken DRGs and quantitatively measure neurite outgrowth. The protocol covers embryo preparation, vertebral column exposure, bilateral DRG microdissection, explant culturing, longitudinal phase-contrast imaging, and neurite image analysis. Results: DRG explants maintained their overall ganglion morphology and remained attached to the culture surface throughout the five-day culture period. Untreated control explants exhibited limited neurite outgrowth, whereas nerve growth factor (NGF) treated explants developed progressively longer and denser radial neurite networks. Neurite outgrowth was readily visualized by live phase contrast microscopy and quantified using an Image-based neurite tracing method, demonstrating the responsiveness of the explant model to NGF and its suitability for quantitative assessment of neurite outgrowth. Conclusion: This simple, low-cost protocol provides an in vitro assay platform for neurotoxicity, neuroregeneration, drug screening, and biomaterial evaluation.

Reported relative biological effectiveness (RBE) values for low-energy X-rays disagree, assays scoring initial DNA double-strand breaks (DSBs) returning about 1.1 and chromosome-level assays 2 to 4. Whether radiation quality varies within the diagnostic range, and how its comparison with a megavoltage reference depends on target scale, has not been quantified on a tube-potential series. A tungsten-anode tube with 1 mm Be and 2.5 mm Al filtration, with copper added in some cases, was modelled in PHITS for 40 to 200 kV. The spectra were transported into a water phantom in which absorbed dose, lineal-energy densities and cluster size distributions were scored for target diameters of 3 nm to 1 micrometre against a cobalt-60 reference; DSB yields were computed in the electron track-structure mode with the PHITS DNA damage tally. Between 40 and 120 kV the depth-dose ratio changed by a factor of 5.7 and the tube output by a factor of 42, whereas the dose-mean lineal energy varied by 2.5 % at 1 micrometre and 1.2 % at 3 nm against a reproducibility of 0.3 %. Relative to cobalt-60 it was 2.05 times larger at 1 micrometre but only 1.08 times larger at 3 nm, while DSB yields per unit dose were 5 to 7 % higher and constant across the range within the 2 % bound set by the statistics. Tube potential therefore changes the amount and distribution of dose but not its physical quality, and a stated RBE is incomplete without the target scale implied by the endpoint.

The chick embryo is a widely used vertebrate model for developmental biology, cardiovascular, and bioengineering research. We present an automated laser speckle contrast imaging (LSCI) platform for non-invasive, longitudinal, label-free imaging of the developing heart and extraembryonic vasculature of the chick embryo through the intact eggshell. The platform can detect blood vessels as small as 70 +/- 17 um in diameter and enable repeated imaging from day 0 to day 8 of incubation without disturbing normal embryonic development. It provides both structural (a map of growing extraembryonic blood vessels) and functional (blood flow dynamics and heart rate) information. Using this system, we generated the first continuous, non-invasive, time-lapse recordings of the early functioning chick heart and its extraembryonic vessels. Longitudinal imaging was performed on 15 embryos from different chicken breeds with eggshell colors white, brown and green. Analysis revealed slight variability in early cardiac and vascular morphogenesis across breeds and shell types and shows its potential for automated developmental staging. Finally, for the first time, the growing vasculature of the yolk sac and the chorioallantoic (CAM) membrane was imaged from day 4 to 8. The heart rate and growth of the extraembryonic blood vessels were quantified automatically during imaging.

Behavioral test batteries are increasingly used to characterize multiple functional domains in zebrafish, yet the potential impact of test sequence on behavioral outcomes remains poorly defined. Here, we systematically evaluated whether test order influences behavioral responses in a three-assay battery comprising the novel tank test (NTT), mirror-induced aggression (MIA), and social preference (SP) test. Adult zebrafish (Danio rerio) were exposed to all possible permutations of the three assays in a fully counterbalanced design, allowing assessment of order effects across locomotor, anxiety-like, aggression-related, and social behaviors. Test order produced modest and parameter-specific effects, primarily affecting locomotor activity in the NTT and social proximity in the SP assay. Time-course analysis revealed within-test behavioral dynamics, with limited evidence that test order modulates early adaptation or late engagement with the testing environment but does not alter overall temporal response profiles. Sex-dependent effects were assay-specific and most pronounced in the NTT, with no consistent sex differences observed in MIA or SP. To evaluate the global structure of behavioral variation, Principal Component Analysis (PCA) was performed across assays. Despite localized effects of test order, no clear multivariate separation between test sequences was observed, indicating that sequential testing does not produce distinct baseline phenotypes. Together, these findings support the robustness and reproducibility of multidomain behavioral batteries while highlighting the importance of standardized test-order reporting to improve cross-study comparability.

The neural system at the periphery is a frontline for embodied cognition, yet an acute, mild perturbation to dissect functional causality is hard to achieve. Peripheral neural processes and the central nervous system may work in concert to generate sophisticated representations regarding self and environments in the brain. This hypothesis, together with the underlying mechanisms, is particularly difficult to test for certain sensory inputs due to the lack of reversible manipulation techniques. Long postulated as a component for path integration, proprioception is one of such modalities. In this study, we developed a murine experimental system to manipulate proprioceptive inputs during memory-dependent localization task (which required precise operant-conditioned licks) in spatial virtual reality (VR). Through bioluminescent optogenetics (luminopsin) selectively expressed in the parvalbumin-positive neurons of the dorsal root ganglia in mice, proprioceptive processing was compromised directly from the periphery to bypass the bottleneck of specific central targeting, which results from the lack of anatomically or genetically dedicated proprioceptive circuits in the brain. In-vivo IVIS imaging and behavior suggested the effects of luminopsin last for roughly 20 minutes. While mice exhibited normal performance in tasks relying on gross motor skills, they showed subtle deficits in challenging spatial tasks that required integration of past movements. These observations support a task-specific role for proprioception, and demonstrate a potential of chemogenetics-like, rapidly reversible strategies for characterizing peripherally defined sensory contribution to spatial cognition. Future work will optimize this approach; for instance, to activate opsins by light with millisecond precision. To our knowledge, this is a first causal demonstration for acute participation of proprioception in path-integration biomechanics, enabling the first temporally defined method for mild perturbation of path-integration mechanisms.

Body odors inform human kin recognition. However, underlying chemical mechanisms and how they relate to body odor perception are still unclear. The present study therefore simultaneously examined perceptual body odor similarity according to human raters, its chemical basis and the link between the two in human families. The study included 33 raters who assessed body odors collected from 60 odor donors belonging to 15 biological families. Each family comprised a mother, a father and one adolescent same-sex sibling-dyad (8 sister dyads, 7 brother dyads). Body odor sampling was conducted with cotton pads for presentation to human raters and with thermal desorption (TD) tubes for Gas Chromatography-Mass Spectrometry (GC-MS) analysis. Relatedness and sex predicted perceived body odor similarity in parent-offspring versus unrelated control dyads. Analyses of chemical profiles revealed greater body odor similarity within families than between unrelated individuals. Parallel to the perceptual results, relatedness and sex predicted chemical similarity based on whole odor profiles. Finally, a machine-learning prediction model integrating perceptual and chemical data, revealed a small but significant positive association between chemical similarity of body odor and perceived olfactory similarity according to human raters, suggesting that odor perception is partly informed by chemical characteristics. Overall, the present study extends our understanding of human chemosensory communication and supports the notion of a family-specific olfactory signature.

The specificity protein 7 (SP7), also known as Osterix (OSX), is a zinc-finger transcription factor essential for osteoblast differentiation and skeletal development. Hereafter, the protein is referred to as OSX/SP7 throughout the manuscript. Although OSX/SP7 has been extensively studied in mammals and teleost fish, its developmental expression pattern in chicken (Gallus gallus) has not been described. Using an experimentally isolated chicken OSX/SP7 sequence, we examined OSX/SP7 mRNA expression during embryogenesis by in situ hybridization. As expected, OSX/SP7 expression was detected in developing skeletal elements undergoing ossification. Unexpectedly, transcripts were also observed in the neuroepithelium, retina, central nervous system, embryonic muscles and integument. Notably, OSX/SP7 expression was present in the neural tube from Hamburger and Hamilton stage 9 and persisted in the developing central nervous system until at least HH40, suggesting that OSX/SP7 functions during avian development may extend beyond osteogenesis. In addition, we reconstructed the chicken OSX/SP7 coding sequence and inferred its associated untranslated regions. The reconstructed ORF was independently supported by maternal and paternal haplotype-resolved chicken genome assemblies and retained the characteristic domain architecture of vertebrate OSX/SP7 proteins despite substantial divergence outside the DNA-binding domain. Taken together, these findings broaden the developmental landscape of OSX/SP7 expression in birds and provide new molecular resources for future studies of OSX/SP7 regulation and function during vertebrate development.

Gelatinous zooplankton are abundant predators and prey in marine food webs, yet their role in helminth (parasitic worms) transmission remains poorly resolved. Here we combine a global synthesis of published records with new morphological and molecular observations from the Mediterranean, Red, Celtic, Baltic and North Seas to assess the ecological role of helminths associated with gelatinous zooplankton. We compiled 431 host-parasite association records from 89 sources, including 23 new records, and show that helminth occurrence and richness in gelatinous hosts are concentrated at temperate latitudes, contrasting with the classical latitudinal diversity gradient. Our sampling revealed markedly higher parasite prevalence, abundance and diversity in the Red Sea than in the Mediterranean, while no helminths were detected in gelatinous zooplankton sampled from the Baltic and North Seas. All helminths recorded were larval stages, indicating that gelatinous zooplankton function as key intermediate hosts in marine helminth life cycles. Molecular analyses identified cestodes, nematodes and digenean trematodes associated with cnidarians, ctenophores, and chaetognaths, including the first record of trematode larvae in a pelagic tunicate. Our findings challenge the view of gelatinous zooplankton as dead-end hosts and identify them as overlooked vectors that may shape marine parasite biogeography, food web connectivity and invasion dynamics.

Background: Determining the cause of death in electrocution is often challenging at autopsy, particularly when characteristic external electrical marks are absent. Reliable internal markers are scarce, and the neuropathological changes associated with electrocution remain underexplored. Objective: To identify and characterize histopathological alterations in the brain in cases of fatal electrocution and evaluate their diagnostic utility. Methods: A prospective observational study was conducted on 12 autopsy cases, including 7 confirmed electrocution deaths and 5 non-electrocution deaths (control). Brain tissue samples from multiple regions were processed using H&E and Golgi Cox staining. Histopathological parameters, including neuronal morphology, vascular changes, haemorrhage, and dendritic architecture, were evaluated. Results: Electrocution cases demonstrated neuronal destruction, axonal ballooning, focal hemorrhages, and dendritic disruption. Cerebellar involvement with Purkinje cell damage. Conclusion: Comprehensive neurohistological technique may helpful in evaluation neuronal damage in electrocution death. Neuropathological findings, particularly neuronal and dendritic injury, may serve as supportive evidence in determining cause of death when external signs are absent. Keywords: Electrocution, Neuropathology, Golgi Cox stain, Forensic histopathology, death.

Electrical stunning is a humane euthanasia method for larval zebrafish, and a suitable euthanasia protocol has recently been established. Here, we assess the associated stress response and tissue integrity. When compared to MS-222 overdose and hypothermia, electrical stunning caused no significant Hsp70 increase. Brain tissue was preserved and muscle tissue affected only slightly. Our results support the suitability of electrical stunning for humane euthanasia and artifact-free post-mortem analyses.

Numerous studies have been conducted to investigate touch on the hand and forearm, but less is known about the rest of the body, particularly the torso and breasts, which play an important role in affective touch. We conducted psychophysical tests on three skin sites on the torso: upper breast, nipple-areola complex (NAC), and under the breast, as well as on two control sites: palm and forearm. We studied touch pleasantness and intensity by applying gentle caress-like strokes using a soft brush, and measured temperature and tactile detection thresholds. The results show an inverted U-shaped relationship for pleasantness over increasing stroking velocity for all skin sites, although we found some differences between sites, particularly at 30 cm.s-1, where participants found stroking on the NAC to be more intense and women additionally found it to be less pleasant. For temperature detection, we found that the NAC was more sensitive in men. For tactile detection, we mainly observed that the areola was far less sensitive than the other areas. Together, these results provide insights into tactile perception in the torso on a fundamental level, which can be applied to understanding how to regain sensations after mastectomy.

Radiotherapy is common practice to treat cancer but produces significant side effects such as chronic pain. Cancer survivors report developing chronic pain due to their treatment even long after the cancer is cured. To understand the mechanisms underlying the radiotherapy-induced chronic pain, we assessed how ionising X-ray radiation exposure during 4 consecutive days of 5 Gy (total radiation dose of 20 Gy) affected dorsal root ganglia (DRG) sensory neurons (rodent F11 cell line). On the 5th day, we assessed known impacts of ionising radiation (senescence, oxidative stress, cellular metabolism, mitochondrial copy number, and mitochondrial respiration) followed by assessing expression of genes associated with populations of DRG neuronal fibres. We discovered that fractionated exposure to ionising radiation increased senescence, mitochondrial copy number, and modulated the NAD+/NADH pathway, but did not change the oxygen consumption rate nor induce oxidative stress 24 hours after the last irradiation exposure. Additionally, ionising radiation altered the expression of genes associated with mechanoreceptor fibres, known to have pro-nociceptive properties in the context of injury and chronic pain.

Despite the well-known health risks of neutron exposures, key gaps remain in understanding neutron-induced molecular responses and identifying reliable biodosimetric markers that distinguish neutrons from photon exposure. We provide the first genome-wide analysis of the human blood transcriptional response to an accelerator-derived fission-like spectrum of neutrons versus photons, evaluating transcriptomic relative biological effectiveness (RBE) and radiation quality-discriminating gene signatures. Whole blood from healthy donors was irradiated ex vivo with X-rays (140 kV, 0-4 Gy, n = 3) or neutrons (0.1-8 MeV, 0-1 Gy, n = 2), incubated for 6 h or 24 h, and processed for RNA sequencing from peripheral blood mononuclear cells (PBMCs). Neutrons were markedly more potent than X-rays at inducing differentially expressed genes (DEGs) at equal doses, showing a peak response 6 h post-irradiation followed by a decline. In contrast, X-rays caused a continuous increase in DEGs up to 24 h (neutrons vs. X-rays at 1 Gy: 1,449 vs. 121 DEGs at 6 h; 996 vs. 621 DEGs at 24 h). A universal p53-centered 34-gene signature, including FDXR, EDA2R, GADD45A, and ZMAT3, showed highly monotonic dose responses (Spearman correlation coefficient {approx} 1) across donors, radiation qualities, and timepoints. Additionally, difference-in-differences analysis identified radiation quality-discriminating genes only at 6 h, with transcriptional convergence observed by 24 h, suggesting a very narrow time window for biodosimetric differentiation. We identified a neutron-specific gene signature driven by cGAS-STING-NF-{kappa}B signaling (RELB, NFKB1, C3, MALAT1) and suppression of B-cell and myeloid identity genes (IGHD, TCL1A, CLEC7A, TLR2), defining a biologically coherent neutron quality index with distinct immunomodulatory effects. For the first time, we assessed neutron RBEs at the gene, pathway, and global transcriptomic levels in a human blood model, reporting a global transcriptomic neutron RBE of 1.30 (95% CI: 1.14-1.49) at 6 h and 1.21 (95% CI: 1.14-1.28) at 24 h, providing a valuable basis for biodosimetry in mixed-field exposure scenarios. Our findings advance the mechanistic understanding of neutron radiation responses and support the development of biodosimetric approaches for mixed-field exposure scenarios.